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Abstract. Employing the perturbative QCD inspired parton cascade, we investigate 

' kinetic and chemical equilibration of the partonic matter created in central heavy ion 

collisions at RHIC and LHC energies. Two types of initial conditions are chosen. One 

is generated by the model of wounded nucleons using the PYTHIA event generator 

^ ' and Glauber geometry. Another is considered as a color glass condensate. We 

show that kinetic equilibration is almost independent on the chosen initial conditions, 

whereas there is a sensitive dependence for chemical equilibration. The time scale 

^ ' of thermalization lies between 1 and 1.5 fm/c. The final parton transverse energy 

0^ ! obtained from BAMPS calculations is compared with the RHIC data and is estimated 

' for the LHC energy. 
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OO ! 1. Introduction 

o 

. ^ ' Comparisons between tlie calculated elliptic flow f 2 from ideal (viscous) hydrodynamics 
^ ' [HE] and the measured f 2 from the experiments [3] at the BNL Relativistic Heavy Ion 

■ Collider (RHIC) support the creation of a nearly equilibrated quark-gluon system flowing 
with a small shear viscosity to entropy density ratio rj/s. On the other hand, initially 
produced quarks and gluons are far from thermal equilibrium due to the subsequent 
rapid longitudinal expansion. How quarks and gluons thermalize within a short time 
scale < 1 fm/c as assumed in hydrodynamical calculations is an important issue. This 
is not only because the thermalization time scale has to be theoretically determined, 
but also because the mechanism that drives the system toward equilibrium should also 
respond to the smallness of the rj/s value and the buildup of the collective flow of the 
quark gluon plasma (QGP). 

To investigate thermalization and collectivity in a spatially expanding particle 
system such like the situation in ultrarelativistic heavy ion collisions, a 3+1 dimensional 
parton cascade Boltzmann Approach of MultiParton Scatterings (BAMPS) [1] is 
developed. In this talk we will first review the role of perturbative QCD (pQCD) 
gluon bremsstrahlung {gg ^ ggg) in thermal equilibrium and in flow buildup. Second, 
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new results on the initial condition dependence of thermalization and decrease of the 
transverse energy are given. Previous studies can be found in [5],|6]. 

We mention that coherent quantum effects like color instabilities [7j may play a 
role in isotropization of particle degrees of freedom at the very initial stage where the 
matter is super dense. However, more quantitative studies are needed to determine their 
significance on the true thermal equilibration in the expanding quark gluon matter. 

2. Parton Cascade BAMPS 

BAMPS solves the Boltzmann equation for partons with pQCD interactions, which 
include, for the moment, gluon elastic scatterings and gluon bremsstrahlung and its 
backreaction. The structure of BAMPS is based on the stochastic interpretation of the 
transition rate [4j, which ensures full detailed balance for multiple scatterings. The 
recent numerical setup can be found in [8]. The critical energy density is set to Cc = 0.6 
GeVfm-l Gluons are terminated when the local energy density is smaller than Cc- 
The differential cross sections and the effective matrix elements are given by ^ 

~ {ql + mlY ' ^> 
I.. 1 2 95'"' 12g'^q\ . i \/ n 

where g'^ = Anas- q± and denote the perpendicular component of the momentum 
transfer and of the radiated gluon momentum in the center-of-mass frame of the 
collision, respectively, y is the momentum rapidity of the radiated gluon in the center- 
of-mass frame, and is the gluon total mean free path, which is calculated self- 
consistently [1]. The interactions of the massless gluons are screened by a Debye 
mass m|) = ndcasNcJ d^p/ {2n)^ f /p where do = 16 is the gluon degeneracy factor 
for Nc = 3. nif) is calculated locally using the gluon density function / obtained from 
the BAMPS simulation. 

The suppression of the bremsstrahlung due to the Landau-Pomeranchuk-Migdal 
(LPM) effect is taken into account within the Bethe-Heitler regime employing the step 
function in equation ([2]). The mean free path serves here as an infrared regulator, 
which leads to a lower cutoff for the transverse momentum of the radiated (or absorbed) 
gluon. Compared to elastic collisions, the collision angles in a bremsstrahlung become 
larger due to the additional regulator [4]. This makes the gg ^ ggg processes much 
more efficient for kinetic equilibration. Only for large value of A^a/s, in an ultrahigh 
energy collision for instance, the cutoff is small and the distribution of collision angles 
is forwards directed, which is similar to the elastic case jlOj . 

3. The role of gg ^ ggg in thermalization and in flow buildup 

The left panel of Fig. 1 shows the momentum isotropization obtained from BAMPS 
calculations with a constant QCD coupling of as = 0.3 for a central Au+Au collision at 
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■s/snn = 200 GeV. The initial gluon distributions are taken as an ensemble of minijets 
with transverse momenta greater than 1.4 GeV [lOj, produced via semihard nucleon- 
nucleon collisions with the Glauber geometry. The results of {p^/E"^) are extracted at 
the collision center within the space time rapidity —0.2 < rjs < 0.2 and transverse radius 
xt < 1-5 fm. Whereas the only elastic pQCD scatterings cannot drive the system toward 
equilibrium, the inclusion of the gluon bremsstrahlung and its backreaction enormously 
speeds up the kinetic equilibration. 

The large-angle distribution in bremsstrahlung is the reason for its efficiency in 
thermalization. Quantitatively we demonstrated that the contributions of different 
processes to momentum isotropization are quantified by the transport rates [10] 



= 71 TTJT; ' \'^) 



n I 



\-m) 



where Ci[/], functional of the gluon density distribution f{p,x), is the corresponding 
collision term describing various interactions, i = gg ^ 99,99 999,999 ~^ 99, 
respectively. The sum of Rf and that of particle drift gives the inverse of the time scale of 
momentum isotropization [TU], which also marks the time scale of overall thermalization. 
The right panel of Fig. 1 shows the transport collision rate, scaled by temperature 
T = 1//?, for elastic gg gg scattering and bremsstrahlung gg ggg, respectively. 
R%^ggg IS & factor of 3 — 5 larger than R^gg^gg over a range in the coupling constant 
from 10~'^ to 0.8, which demonstrates the essential role of the bremsstrahlung in thermal 
equilibration. 

In addition, the shear viscosity to the entropy density ratio t]/s is inversely 
proportional to the sum of transport collision rate [1^. Thus gg ^ ggg processes 
significantly decrease the rj/s value. We also found [12] that gg ggg processes build 
up large elliptic fiow t>2 observed at RHIC. When ag = 0.3 — 0.6 is chosen, the calculated 
V2 from BAMPS is comparable with the experimental data and the extracted QGP rj/s 
lies between 0.15 for a. = 0.3 and 0.08 for a„ = 0.6. 
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4. Initial conditions: model of wounded nucleons and color glass condensate 

To study the initial condition dependence of thermalization we choose two different 
types of parton initial conditions in nucleus-nucleus collisions. One is based on the 
model of wounded nucleons with the Glauber geometry [4]. A nucleus-nucleus collision 
is considered as sequent binary nucleon-nucleon collisions. To obtain the parton 
momentum distribution in a proton-proton collision we employ the PYTHIA event 
generator [13] and turn down the function of the parton fragmentation. On shell 
quarks and gluons are produced either by (semi)-hard parton-parton collisions or by 
the associated initial and final state radiations. Soft partons produced are not included 
to cascade calculations. To obtain initial conditions in a central Au-|-Au (Pb-|-Pb) 
collision the parton production in a p+p collision is scaled by a number of binary 
collisions A^bin, which is set to be A^bin = 1000. The positions of the initial partons are 
determined according to the nuclear overlapping density within the Glauber geometry 
using the Woods-Saxon profile [ID]. We note that exact fractions of proton-proton, 
proton-neutron, and neutron-neutron collisions and shadowing effects in a nucleus- 
nucleus collision will be taken into account in a forthcoming paper. 

Another type of initial conditions is considered as the color glass condensate (CGC) 
within a /cT^-factorization KLN approach [151 [IB] 

dNg AN^ fPT'' d^pT 



J c?2/i;Tas0A(a;i,k|;rr)0B(x2,(pr-kr)^;rr) 



d^rxdy — 1 J pj, 

Applying such initial conditions to ideal hydrodynamic calculations final hadronic yields 
and their spectra measured at RHIC are well reproduced [U] . In addition, CGC initial 
conditions give larger initial eccentricity than the Glauber-type ones [16]. This leads 
to larger elliptic flow in noncentral nucleus- nucleus collisions [2]. The initial condition 
dependence on V2 from BAMPS calculations will be investigated in the near future. 
In this work, the unintegrated gluon distribution r^) is taken from Ref. |16j . 

which, compared to the original KLN approach, gives a smooth transition from the 
saturation, (f){k'^) = const., towards the DGLAP regime, (f){k^) ~ 1/k'^. We set the 
prefactor for k^; r^) so that the total gluon energy is 80% of y/s = 200(5500)A GeV 
at RHIC(LHC), the same as obtained for quarks and gluons together in the wounded 
nucleons model. 



5. Energy decrease and thermalization at RHIC and LHC: dependence on 
initial conditions 

The left panel of Fig. 2 shows the momentum rapidity distribution of transverse energy 
of initial quarks and gluons produced in wounded nucleons model (wn) and CGC 
approach in a central Au-|-Au collision at RHIC. The transverse energy of CGC gluons 
is larger than that of quarks and gluons from wn model over a wide range of rapidity 
except for at large rapidity \y\ > 4. The difference stems from the different approach 
for the production of gluons with low transverse momentum: Whereas the CGC gives 
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Figure 2. (color online) Left panel: momentum rapidity distribution of transverse 
energy of initial quarks and gluons in a central Au+Au collision at RHIC. Right panel: 
time evolution of the total transverse energy at midrapidity. 



a moderate yield (even though suppressed due to the saturation), no production is 
expected in the wn model below a scale of = 1-4 — 2 GeV, below which the pQCD 
parton scatterings is no longer valid. 

The right panel of Fig. 2 shows the decrease of the transverse energy at midrapidity 
calculated from BAMPS simulations using = 0.3. For both types of initial conditions 
a formation time for each parton is introduced as tq = 0.15 cosh?/ fm/c, where 0.15 fm/c 
is the overlapping time of a Au+Au collision. The solid curve depicts the result when 
quarks in the wn model do not interact, while the dashed curve depicts the result when 
quarks interact as strong as gluons. The implementation of real pQCD quark dynamics 
in BAMPS calculations is in progress. We expect that the true final transverse energy of 
quarks and gluons will lie between the final values of the solid and dashed curve and thus 
will be comparable with the exprimental data including hadronic and electromagnetic 
components, cIEt / dy\y=o = 620 ± 33 GeV [IT] . 

The final transverse energy of CGC gluons is about 30% larger than the 
experimental data due to the larger initial value compared with that in the wn model, 
although the Et decrease of CGC gluons is slighly stronger. To match the experimental 
data, a smaller critical energy density Cc has to be taken to have a later kinetic freezeout, 
or a larger coupling ag has to be assumed, which will lead to a smaller rj/s of the 
QGP. Although the latter will contradict the findings from recent viscous hydrodynamic 
calculations [2], new analyses with experimental data show that the extracted rj/s with 
CGC as the initial conditions is indeed smaller than that with the Glauber-type approach 
|18j . Further comparisons between different extraction models have to be made. 

The local momentum isotropization, {p\/ E'^){t), is shown in the left panel of Fig. 
3. The solid curve depicts the momentum isotropization of gluons in the wn model, 
while the dashed curve depicts the result including quarks that are assumed to interact 
as gluons. Both curves are almost identical. This indicates that the sum of the transport 
collision rates i?**" ([3]), which determines momentum isotropization [10], is the same in 
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Figure 3. (color online) Momentum isotropization (left) and time evolution of the 
fugacity (right) in the central region in central Au+Au collisions at RHIC energy. 

the QCD matter with or without quarks. 

The momentum isotropization of the CGC gluons (dotted curve) is shghtly faster 
than that of partons in the wn modeL Since the initial energy density of the CGC 
gluons is larger (see the left panel of Fig. 2), the local effective temperature T = e/(3n) 
is also larger than that in the wn model. This leads to a larger R^^ , because i?*'" is 
approximately proportional to T [10]. In any case, the momentum isotropization with 
the initial conditions from the wn and CGC approach is not much different from that 
with minijets initial conditions shown in the left panel of Fig. 1. The time scale of the 
momentum isotropization (kinetic equilibration) is about 1.5 fm/c. 

The chemical equilibration is described by the time evolution of the fugacity 
\{t) = n(t)/neq(t), where n is the local particle density and Ueq = dcT^/'n''^ is the 
value at thermal equilibrium. The chemical equilibration shown in the right panel of 
Fig. 3 is quite different between the wn and CGC approach. Whereas quarks and gluons 
in the wn model is initially almost in chemical equilibrium, the CGC gluons need a time 
of 2 fm/c to achieve the chemical equilibrium. 

The rapidity distribution of the initial transverse energy in central Pb+Pb collisions 
at Large Hadron Collider (LHC) is shown in the left panel of Fig. 4. The difference 
between two types of initial conditions is smaller compared to the case at RHIC, because 
the saturation scale Qg of the CGC becomes larger at higher energy collisions and gluon 
production below Qs is suppressed. 

We perform BAMPS calculations using = 0.2 and the formation time tq = 
0.15 cosh ?/ fm/c. The final transverse energy at midrapidity is estimated to be between 
1620 and 2150 GeV (see the right panel of Fig. 4). The notion "wn" in the figure implies 
calculations including quarks, which are assumed to behave as gluons. Note that the 
choices for and tq are crucial assumptions. With larger ag and smaller tq the decrease 
of the transverse energy becomes stronger. 

Figure 5 shows the kinetic (left panel) and chemical (right panel) equilibration at 
the collision center of central Pb+Pb collisions at LHC. Similar to the case at RHIC, 
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Figure 4. (color online) Same as Fig. 2 but for central Pb+Pb collisions at the LHC 
energy ^ = 5500A GeV with = 0.2. 



0.40 
0.35 
0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 







wn 




cgc 




Pb+Pb@LHC 







3 4 5 

t (fm/c) 




2 3 4 5 

t (fm/c) 



Figure 5. (color online) Same as Fig. 3 but for central Pb+Pb collisions at the LHC 
energy with = 0.2. 



there is no much difference in momentum isotropization, whereas the CGC gluons reach 
the chemical equilibrium at a later time. Overall thermal equilibrium is expected to be 
achieved at about 1 fm/c according to the relaxation formula for {p\/ E'^) |lUj . 

6. Summary 

Employing a parton cascade BAMPS we investigated thermalization of the QCD matter 
produced in central nucleus-nucleus collisions at the RHIC and LHC energy. The 
transport collision rate is a proper quantity describing kinetic equilibration and separates 
contributions from various interactions. We found the QCD inspired bremsstrahlung is 
the dominant process in thermal equilibration and in flow buildup. 

We studied the dependence of thermalization on initial conditions within the 
wounded nucleons model and the CGC approach. Using the same value of the 
momentum isotropization is almost independent on the initial conditions, whereas CGC 
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gluons reach chemical equihbrium at a later time than partons in the wn model. On 
the other hand, the final transverse energy of CGC gluons is higher than that in the 
wn model due to their initial difference. A larger as for CGC gluons will accelerate the 
decrease of the transverse energy to reach a comparable value with the experimental 
data. This will also accelerate thermalization of CGC and will lead to a smaller r]/s 
ratio. Detailed study including real quark dynamics and connecting the initial condition 
dependence of the buildup of elliptic fiow f 2 will be done in the near future. 
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